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For a regular sub-Riemannian manifold we study the Radon-Nikodym derivative of the spher- 
ical HausdorfF measure with respect to a smooth volume. We prove that this is the volume of 
the unit ball in the nilpotent approximation and it is always a continuous function. We then 
prove that up to dimension 4 it is smooth, while starting from dimension 5, in corank 1 case, it 



q 

(~| , is C^ (and C* on every smooth curve) but in general not C^. These results answer to a question 



addressed by Montgomery about the relation between two intrinsic volumes that can be defined 
in a sub-Riemannian manifold, namely the Popp and the HausdorfF volume. If the nilpotent 
approximation depends on the point (that may happen starting from dimension 5), then they 
are not proportional, in general. 

cn 
>, 
Q ■ 1 Introduction 

in 

f"^ ■ In this paper, by a sub-Riemannian manifold we mean a triple S = (M, A,g), where M is a 

iy-> , connected orientable smooth manifold of dimension n, A is a smooth vector distribution of constant 

^^ ' rank k < n, satisfying the Hormander condition and g is an Euclidean structure on A. 



A sub-Riemannian manifold has a natural structure of metric space, where the distance is the 
so called Carnot-Caratheodory distance 



rS ; d{qo, qi) = inf{ / yg^(t){iit),'yit)) dt \ -/: [0,T] ^ M is a Lipschitz curve, 

^ ' ° 7(0) = qo,7{T) = qu i{t) G A^(j) a.e. in [0,r]}. 



As a consequence of the Hormander condition this distance is always finite and continuous, and 
induces on M the original topology (see ||2l|, p9[). 

Since (M, d) is a metric space, for every a > one can define the a-dimensional Hausdorff 
measure on M, and compute the Hausdorff dimension of M. 

Define A^ := A,A*+^ := A* -|- [A*, A], for every i = 1,2,... Under the hypothesis that the 
sub-Riemannian manifold is regular, i.e. if the dimension of A', i = 1, . . . , ?n, do not depend on 
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the point, the Hormander condition guarantees that there exists (a miminial) m G N, called step 
of the structure, such that A^ = TqM, for all q G M. The sequence 

g{S) := (dim A, dim A^,..., dim A"') 

II II 

k " 

is called growth vector of the structure. 

In this case, the graded vector space associated to the filtration Ag C A^ C . . . C A^ = TqM, 



gr,(A) = 0A^/A^-S where A° 



is well defined. Moreover, it is well known that the Hausdorff dimension of M is given by the 
formula (see [p5| ) 



Y,^ki, ki:= dim Ai/Al-\ 



Q = ^iki, /ci :=dimA*g/A: 



In particular the Hausdorff dimension is always bigger than the topological dimension of M. 

Moreover, the Q-dimensional Hausdorff measure (denoted by 'H'°^ in the following) behaves like 
a volume. More precisely, in [35] Mitchell proved that if /i is a smooth volumq^ on M , then 



dji = ffj.'H.d'H^ , where /^-^ is a positive measurable function that is locally bounded and locally 
bounded away from zero, that is the Radon-Nikodym derivative of fj, with respect to 'H'^. According 
to Mitchell terminology, this means that the two measures are commensurable one with respect to 
the other. 

Notice that the Hausdorff measure on sub-Riemannian manifolds has been intensively studied, 
see for instance p^, E5| . A deep study of the Hausdorff measure for hypersurfaces in sub-Riemannian 



geometry, in particular in the context of Carnot groups, can be found in |^, [l^, |ll|, |lj, |l^, |2^, |3^, ^| 



and references therein. Hausdorff measures for curves in sub-Riemannian manifolds were also 



studied in the problem of motion planning and complexity, see [24, ^, 26, ^]. 

Let us recall that there are two common non-equivalent definitions of Hausdorff measure. The 
standard Hausdorff measure, where arbitrary coverings can be used, and the spherical Hausdorff 
measure, where only ball-coverings appear (see Definition p!9| ). However it is well known that, 
if S'^ denotes the spherical Q-Hausdorff measure, then T-L'^ is commensurable with 5*^.F| As a 
consequence, 5*^ is commensurable with ^u, i.e. 

dfi = UsdS"^, 

for a positive measurable function f^s that is locally bounded and locally bounded away from zero. 
In this paper, we are interested to the properties of the function f^s. In particular, we would like 
to get informations about its regularity. 

The reason why we study the spherical Hausdorff measure and not the standard Hausdorff mea- 
sure is that the first one appears to be more natural. Indeed, as explained later, /^^ is determined by 



^In the following by a smooth volume on M we mean a measure /j, associated to a smooth non-vanishing n-form 
ojp G A"M, i.e. for every measurable subset A G M we set /i(A) = /^t^^/j- 

^Indeed they are absolutely continuous one with respect to the other. In particular, for every a > 0, we have 
2-Q,^a < -^o^ < 5°^ (see for instance ii). 



the volume of the unit sub-Riemannian ball of the nilpotent approximation of the sub-Riemannian 
manifold, that can be explicitly described in a certain number of cases (see Theorem |l| below). On 
the other hand nothing is known on how to compute ffj^-u- We conjecture that /^-^ is given by the 
//-volume of certain isodiametric sets, i.e. the maximum of the /i- volume among all sets of diameter 
1 in the nilpotent approximation (see [^, ^] and reference therein for a discussion on isodiametric 
sets). This quantity is not very natural in sub-Riemannian geometry and is extremely difficult to 
compute. 

Our interests in studying /^^ comes from the following question: 

Ql How can we define an intrinsic volume in a sub-Riemannian manifold? 

Here by intrinsic we mean a volume which depends neither on the choice of the coordinate system, 
nor on the choice of the orthonormal frame, but only on the sub-Riemannian structure. 



This question was first pointed out by Brockett, see [0], and by Montgomery in his book |36|. 
Having a volume that depends only on the geometric structure is interesting by itself, however, it 
is also necessary to define intrinsically a Laplacian in a sub-Riemannian manifold. We recall that 
the Laplacian is defined as the divergence of the gradient and the definition of the divergence needs 
a volume since it measures how much the fiow of a vector field increases or decreases the volume. 

Before talking about the question Ql in sub-Riemannian geometry, let us briefiy discuss it in 
the Riemannian case. In a n-dimensional Riemannian manifold there are three common ways of 
defining an invariant volume. The first is defined through the Riemannian structure and it is the 
so called Riemannian volume, which in coordinates has the expression ^dx^ . . . dx^, where g is 
the determinant of the metric. The second and the third ones are defined via the Riemannian 
distance and are the n-dimensional Hausdorff measure and the n-dimensional spherical Hausdorff 
measure. These three volumes are indeed proportional (the constant of proportionality depending 
on the normalization, see e.g. [^, p^]). 

For what concern sub-Riemannian geometry, a regular sub-Riemannian manifold is a metric 
space, hence it is possible to define the Hausdorff volume 'H'^ and the spherical Hausdorff volume 
S^ . Also, there is an equivalent of the Riemannian volume, the so called Popp's volume V, 
introduced by Montgomery in his book ||3^ (see also [^]). The Popp volume is a smooth volume 
and was used in Q to define intrinsically the Laplacian (indeed a sub-Laplacian) in sub-Riemannian 
geometry. 

In his book, Montgomery proposed to study whether these invariant volumes are proportional 
as it occours in Riemannian geometry. More precisely, he addressed the following question: 

Q2 Is Popp's measure equal to a constant multiple (perhaps depending on the growth vector) of 
the Hausdorff measure? 

Mongomery noted that the answer to this question is positive for left-invariant sub-Riemannian 
structures on Lie groups, since the Hausdorff (both the standard and spherical one) and the Popp 
volumes are left-invariant and hence proportional to the left Haar measure. But this question is 
nontrivial when there is no group structure. 

One of the main purpose of our analysis is to answer to question Q2 for the spherical Hausdorff 
measure, i.e. to the question if the function f-ps (defined by dV = f-psdS^) is constant or not. 
More precisely, we get a positive answer for regular sub-Riemannian manifolds of dimension 3 and 
4, while a negative answer starting from dimension 5, in general. 

Once a negative answer to Q2 is given, it is natural to ask 



Q3 What is the regularity of /-p^? 

This question is important since the definition of an intrinsic Laplacian via S^ require /■ps to be 
at least C^. 

Notice that since the Popp measure is a smooth volume, then /^^ is C'^, A; = 0, 1, . . . , cxo if and 
only if f-ps is as well. 

We prove that /^5 is a continuous function and that for n < 4 it is smooth. In dimension 5 it 
is C^ but not smooth, in general. Moreover, we prove that the same result holds in all corank 1 
cases (see Section |^ for a precise definition) . 

Our main tool is the nilpotent approximation (or the symbol) of the sub-Riemannian structure. 
Recall that, under the regularity hypothesis, the sub-Riemannian structure S = (M, A, g) induces 
a structure of nilpotent Lie algebra on gr (A). The nilpotent approximation at q is the nilpotent 
simply connected Lie group Grq(A) generated by this Lie algebra, endowed with a suitable left- 
invariant sub-Riemannian structure Sg induced by S, as explained in Section ^. 

Recall that there exists a canonical isomorphism of 1-dimensional vector spaces (see Q for 
details) 

n n 

/\(r;M):^/\(gr,(A)*). (1) 

Given a smooth volume /i on M, we define the induced volume /Ig on the nilpotent approximation 
at point q as the left-invariant volume on Grg(A) canonically associated to w^(q') E /\^{T*M) by 
the above isomorphism. 

The first result concerns an explicit formula for /^5. 

Theorem 1. Let S = (Af, A,g) he a regular sub-Riemannian manifold. Let fi a volume on M and 
Jiq the induced volume on the nilpotent approximation at point q £ AI . If A is an open subset of 
M, then 

KA) = ^lj,{B,)dS'^, (2) 

where Bq is the unit hall in the nilpotent approximation at point q, i. e. 

Starting from this formula we prove our first result about regularity of the density: 

Corollary 2. Let S = (M, A,g) be a regular sub-Riemannian manifold and let ^ be a smooth 
volume on M . Then the density f^s is a continuous function. 

Theorem |l|, specified for the Popp measure V, permits to answer the Montgomery's question. 
Indeed, the measure Vq induced by V on the nilpotent approximation at point q coincides with the 
Popp measure built on S^, as a sub-Riemannian structure. In other words, if we denote V^ the 
Popp measure on Sg, we get 

r, = Vq. (3) 

Hence, if the nilpotent approximation does not depend on the point, then f-ps is constant. In 
other words we have the following corollary. 



Corollary 3. Let S = (M, A, g) be a regular sub-Riemannian manifold and Sq its nilpotent ap- 
proximation at point q G M. // S^^ is isometric to 8^2 for any qi,q2 G M, then f-ps is constant. 
In particular this happens if the sub-Riemannian structure is free. 

For the definition of free structure see p^] . 

Notice that, in the Riemannian case, nilpotent approximations at different points are isometric, 
hence the Hausdorff measure is proportional to the Riemannian volume (see pO, p^ ) . 

When the nilpotent approximation contains parameters that are function of the point, then, in 
general, f-ps is not constant. We have analyzed in details all growth vectors in dimension less or 
equal than 5: 

- dimension 3: (2,3), 

- dimension 4: (2,3,4), (3,4), 

- dimension 5: (2,3,5), (3,5), (4,5) and the non generic cases (2,3,4,5), (3,4,5). 

In all cases the nilpotent approximation is unique, except for the (4,5) case. As a consequence, we 
get: 

Theorem 4. Let S = (Af, A,g) be a regular sub-Riemannian manifold of dimension n < 5. Let fi 
be a smooth volume on M and V be the Popp measure. Then 

(i) ifQiS) 7^ (4,5), then f-pg is constant. As a consequence f^g is smooth. 

(a) ifG{S) = (4,5), then /^^ is C^ (and C^ on smooth curves) but not C^ , in general. 

Actually the regularity result obtained in the (4,5) case holds for all corank 1 structures, as 
specified by the following theorem. 

Theorem 5. Let S = (M, A, g) be a regular corank 1 sub-Riemannian manifold of dimension 
n> 5. Let fi be a smooth volume on M . Then f^g is C^ (and C^ on smooth curves) but not C^ , in 
general. 

Recall that for a corank 1 structure one has ^(S) = (n — l,n) (see also Section 0). 

Notice that Theorem |5| apply in particular for the Popp measure. The loss of regularity of 
ffj^s is due to the presence of what are called resonance points. More precisely, the parameters 
appearing in the nilpotent approximation are the eigenvalues of a certain skew-symmetric matrix 
which depends on the point. Resonances are the points in which these eigenvalues are inultiple. 

To prove Theorem g, we have computed explicitly the optimal synthesis (i.e. all curves that 
minimize distance starting from one point) of the nilpotent approximation and, as a consequence, 
the volume of nilpotent balls Bg. 

Another byproduct of our analysis is 

Proposition 6. Under the hypothesis of Theorem ||, if there are no resonance points then /^5 is 
smooth. 

The structure of the paper is the following. In Section || we recall basic facts about sub- 
Riemannian geometry and about Hausdorff measures. In Section [^ we provide normal forms for 
nilpotent structures in dimension less or equal than 5. In Section ^ we prove Theorem || and its 
corollaries, while in Section |5| we study the differentiability of the density for the corank 1 case. In 
the last Section we prove Theorem ^ 



2 Basic Definitions 

2.1 Sub-Riemannian manifolds 

We start recalling the definition of sub-Riemannian manifold. 

Definition 7. A sub-Riemannian manifold is a triple S = (M, A,g), wfiere 

(i) M is a connected orientable smootfi manifold of dimension n > 3; 

(ii) A is a smooth distribution of constant rank k < n satisfying the Hormander condition, i.e. a 
smooth map that associates to q £ M a /c-dimensional subspace Ag of TqM and we have 

span{[Xi, [. . . [X,„i, Xj]]]{q) | X, G A, j G N} = T,M, ^ q e M, (4) 

where A denotes the set of horizontal smooth vector fields on M, i.e. 

A = {X e Vec(M) I X{q) G Ag V g G M} . 

(Hi) gg is a Riemannian metric on Ag which is smooth as function of q. We denote the norm of a 
vector i; G Ag with |f|, i.e. \v\ = ^/gqivjv). 

A Lipschitz continuous curve 7 : [0, T] — t- M is said to be horizontal (or admissible) if 

7(t) G A^(j) for a.e. t G [0,r]. 

Given an horizontal curve 7 : [0,T] — )• M, the length of -y is 

Z(7)= rm)\dt. (5) 

The distance induced by the sub-Riemannian structure on M is the function 

d{qo,qi) = inf{/(7) | 7(0) = q'o,7(^) = 91,7 horizontal}. (6) 

The hypothesis of connectedness of M and the Hormander condition guarantees the finiteness and 
the continuity of (i(-,-) with respect to the topology of M (Chow-Rashevsky theorem, see, for 
instance, [^]). The function d(-,-) is called the Carnot-Caratheodory distance and gives to M the 
structure of metric space (see |l^, ^). 

Remark 8. It is a standard fact that /(7) is invariant under reparameterization of the curve 7. 
Moreover, if an admissible curve 7 minimizes the so-called action functional 

J{l):=\f^m)\'dt. 

with T fixed (and fixed initial and final point), then |7(t)| is constant and 7 is also a minimizer of 
/(•). On the other side, a minimizer 7 of /(•) such that |7(i)| is constant is a minimizer of J(-) with 
T = l{j)/v. 



Locally, the pair (A, g) can be given by assigning a set of k smooth vector fields spanning A 
and that are orthonormal for g, i.e. 

Ag=span{Xi{q),...,Xkiq)}, Sq{MQ),Xjiq)) = S^j. (7) 

In this case, the set {Xi,...,Xfc} is called a local orthonormal frame for the sub-Riemannian 
structure. 

Definition 9. Let A be a distribution. Its flag is the sequence of distributions A^ C A^ C . . . 
defined through the recursive formula 

A^:=A, A^+^ := A^ + [ASA]. 

A sub-Riemannian manifold is said to be regular if for each i = 1,2,... the dimension of A* 
does not depend on the point qq £ M. 

Remark 10. In this paper we always deal with regular sub-Riemannian manifolds. In this case 
Hormander condition can be rewritten as follows: 

3 minimal m G N such that A™ = TqM, Mq^M. 

The sequence ^(S) := (dim A, dim A^, . . . jdimA*") is called growth vector. Under the regularity 
assumption ^(S) does not depend on the point and m is said the step of the structure. The minimal 
growth is (A;, A; + 1, /c -|- 2, . . . , n). When the growth is maximal the sub-Riemannian structure is 



called free (see |36|). 

A sub-Riemannian manifold is said to be corank 1 if its growth vector satisfies ^(S) = (n — 1, n). 
A sub-Riemannian manifold S of odd dimension is said to be contact if A = ker cj, where ui £ A^M 
and dio\/\ is non degenerate. A sub-Riemannian manifold M of even dimension is said to be quasi- 
contact if A = kerw, where u G A^M and satisfies dimkerdw|A = 1- 

Notice that contact and quasi-contact structures are regular and corank 1. 

A sub-Riemannian manifold is said to be nilpotent if there exists an orthonormal frame for the 
structure {Xi, . . . ,Xfc} and j £ N such that [Xjj, [Xi^, . . . , [Xj^_-^,Xj^.]]] = for every commutator 
of length j. 

2.2 Geodesies 

In this section we briefly recall some facts about sub-Riemannian geodesies. In particular, we define 
the sub-Riemannian Hamiltonian. 

Definition 11. A geodesic for a sub-Riemannian manifold S = (M, A, g) is a curve 7 : [0, T] — )■ M 
such that for every sufficiently small interval [^1,^2] C [0,T], the restriction 71 is a minimizer of 

J(-). A geodesic for which g'y(i)(7(t),7(t)) is (constantly) equal to one is said to be parameterized 
by arclength. 

Let us consider the cotangent bundle T*M with the canonical projection vr : T*M — ?• M, 
and denote the standard pairing between vectors and covectors with (•,•). The Liouville 1-form 
s E A^(T*M) is defined as follows: sx = A o vr^,, for every A E T*M. The canonical symplectic 
structure on T*M is defined by the closed 2-form a = ds. In canonical coordinates (^,x) 



s = ^ Cidxi, a = '^dS,iA dxi. 



i=l i=l 



We denote the Hamiltonian vector field associated to a function h G C°°{T*M) with h. Namely 
we have dh = a{-, h) and in coordinates we have 

-> Y^ dh d dh d 

^ dii dxi dxi d^i 

The sub-Riemannian structure defines an Euclidean norm | • | on the distribution Ag C TqM . 
As a matter of fact this induces a dual norm 



|A|| = max (A, v) , A G T*M, 

D6Aq 



which is well defined on A* ~ T*M/A^, where A^ = {A G T*M\ {X,v) = 0,Vv G Ag} is the 
annichilator of the distribution. 

The sub-Riemannian Hamiltonian is the smooth function on T*M, which is quadratic on fibers, 

defined by 

1 

T 
If {^1, . . . , Xfc} is a local orthonormal frame for the sub-Riemannian structure it is easy to see that 



H{X) = -\\Xf, XeT*M. 



k 
H{X) = -^{X,X,{q)f, XgT*M, q = 7r{X). 



i=l 

Let S = (M, A, g) be a sub-Riemannian manifold and fix qq G M. We define the endpoint map 
(at time 1) as 

F:U^M, F(7) = 7(l), 

where U denotes the set of admissible trajectories starting from qq and defined in [0, 1]. If we fix a 
point qi G M, the problem of finding shortest paths from go to qi is equivalent to the following one 

niin J(7), (8) 

where J is the action functional (see Remark ^). Then Lagrange multipliers rule implies that any 
7 G ZY solution of (P) satisfies one of the following equations 

XiD^F = dyJ, (9) 

XiD^F = 0, (10) 

for some nonzero covector Ai G ^*Q^M associated to 7. The following characterization is a corollary 
of Pontryagin Maximum Principle (PMP for short, see for instance P, fsl, 3C, |38|): 



Theorem 12. Let ^ he a m,inimizer. A nonzero covector Ai satisfies (^) or (|l0|) if and only if 
there exists a Lipschitz curve X{t) G T*,^-.M, t G [0, 1], such that A(l) = Ai and 

- if @j holds, then A(t) is a solution of X{t) = H{X{t)) for a.e. t G [0, 1], 

- if dll) holds, then A(t) satisfies a{X{t),Tx^t)^^) = for a.e. t£ [0, 1]. 



The curve A(t) is said to be an extremal associated to 7(i). In the first case X{t) is called a normal 
extremal while in the second one an abnormal extremal. 

Remark 13. It is possible to give a unified characterization of normal and abnormal extremals in 
terms of the symplectic form. Indeed the Hamiltonian H is always constant on extremals, hence 



A(t) C H~^{c) for some c > 0. Theorem 12 can be rephrased as follows: any extremal X{t) such 
that H(X(t)) = c is a reparametrization of a characteristic curve of the differential form a\ff-ifc\, 
where c = for abnormal extremals, and c > for normal ones. 

Also notice that, if A(t) is a normal extremal, then, for every a > 0, Xa{t) := aX{at) is also 
a normal extremal. If the curve is parametrized in such a way that H = ^ then we say that the 
extremal is arclength parameterized. Trajectories parametrized by arclength corresponds to initial 
covectors Aq belonging to the hypercylinder A^^ := T*^^M nH~^{^) ~ 5^"^ x M"~^ in T*^^M. 

Remark 14. From Theorem |l^ it follows that A(t) = e (Aq) is the normal extremal with initial 
covector Aq G Ag^. If vr : T*M — )• M denotes the canonical projection, then it is well known that 
7(t) = 7r(A(t)) is a geodesic (starting from qq). On the other hand, in every 2-step sub-Riemannian 
manifold all geodesies are projection of normal extremals, since there is no strict abnormal minimizer 
(see Goh conditions, [^). 

The following proposition resumes some basic properties of small sub-Riemannian balls 

Proposition 15. Let S be a sub-Riemannian manifold and Bq^^[£) the sub-Riemannian ball of 
radius e at fixed "point qq £ M . For e > small enough we have: 

(i) Vg G Bqg{e) there exists a minimizer that join q and qq, 

(a) diam(i?g„(e)) = 2e. 

Claim (i) is a consequence of Filippov theorem (see [^, |l^). To prove (ii) it is sufficient to show 
that, for e small enough, there exists two points in qi,q2 S dBqg{£) such that d{qi,q2) = 2e. 

To this purpose, consider the projection 7(t) = 7r(A(t)) of a normal extremal starting from 
7(0) = goi and defined in a small neighborhood of zero t g] — 5, (5[ . Using arguments of Chapter 17 
of [^ one can prove that 7(t) is globally minimizer. Hence if we consider < e < 5 we have that 
qi = 7(— e) and (72 = 7(£) satisfy the property required, which proves claim [ii). 

Definition 16. Fix q^ E M. We define the Exponential map starting from go as 

Expg„ : T;^M ^ M, Expg„(Ao) = 7r(e^(Ao)). 

Using the homogeneity property H{cX) = c'^H{X), Vc > 0, we have that 

e-^(sA) = e*-^(A), Vs>0. 

In other words we can recover the geodesic on the manifold with initial covector Aq as the image 
under Exp^^^ of the ray {tXo,t £ [0, 1]} C T*^M that join the origin to Aq. 

Exp,^(tAo) = vr(e^(tAo)) = ^(e*^(Ao)) = vr(A(t)) = 7(t). 
Next, we recall the definition of cut and conjugate time. 



Definition 17. Let S = (M, A,g) be a sub-Riemannian manifold. Let qq G M and Aq G A^q. 
Assume that the geodesic 7(t) = Expgjj(tAo) for t > 0, is not abnormal. 

(i) The first conjugate time is t(Ao) = min{t > 0, f Aq is a critical point of Exp }. 

{a) The cut time is tc{Xo) = min{t > 0, 3Ai G Ag^, Ai / Aq s.t. Exp^y(tc(Ao)Ao) = Expgjj(tc(Ao)Ai)}. 

It is well known that if a geodesic is not abnormal then it loses optimality either at the cut or 
at the conjugate locus (see for instance Q). 

2.3 Hausdorff measures 

In this section we recall definitions of Hausdorff measure and spherical Hausdorff measure. We 
start with the definition of smooth volume. 

Definition 18. Let M be a n-dimensional smooth manifold, which is connected and orientable. By 
a smooth volum,e on M we mean a measure fj, on M associated to a smooth non-vanishing n-form 
ujf^ £ A"M, i.e. for every subset ^ C M we set 



Ja 



The Popp volume V, which is a smooth volume in the sense of Definition 18, is the volume 



associated to a n-form wp that can be intrinsically defined via the sub-Riemannian structure (see 

Let (M, d) be a metric space and denote with B the set of balls in M. 

Definition 19. Let A be a subset of M and a > 0. 
The a- dimensional Hausdorff measure of A is 

n"iA):=limnnA), 
where 

CX) oo 

nfiA) := inf{^ diam(Ai)", A C \J Ai, diam{Ai) < 6}. 

i=l j=l 

The a- dimensional spherical Hausdorff measure of A is 

S''{A):=^mSf{A), 

where 

oo oo 

S^{A) := inf{^diam(Si)",^ C |J Bi,B, G B,diam{Bi) < 6}. 

These two measures are commensurable since it holds (see |2^) 

2-°cS°(^) <?^"(A) <5°(A), yAcM. (11) 

The Hausdorff dimension of A is defined as 

inf{a > 0, n^iA) = 0} = sup{a > 0, ^"(A) = +oo}. (12) 

10 



Formula ( |ri|) g uarantees that Hausdorff dimension of A does not change if we replace 'H°' with 
5" in formula ([l2|). 

It is a standard fact that the Hausdorff dimension of a Riemannian manifold, considered as a 
metric space, coincides with its topological dimension. On the other side, we have the following 



Theorem 20. Let (il/, A,g) be a regular suh- Riemannian manifold. Its Hausdorff dimension as a 
metric space is 

m 

Q = ^ik„ /ci :=dimA*-dimA*-^ 

Moreover S^ is commensurable to a smooth volume /i on M , i.e. for every compact K C M there 
exists ai, 02 > such that 

aicS'3 <fi< a2S^. (13) 



This theorem was proved by Mitchell in [^5| . In its original version it was stated for the Lebesgue 
measure and the standard Hausdorff measure. 

3 The nilpotent approximation 

In this section we briefly recall the concept of nilpotent approximation. For details see [^, 12|. 



3.1 Privileged coordinates 

Let S = (Af, A,g) be a sub- Riemannian manifold and (Xi,...,Xfc) an orthonormal frame. Fix 
a point q G M and consider the flag of the distribution A^ C A^ C . . . C A™. Recall that 
ki = dim A* — dim A*"-"^ for i = 1, . . . , m, and that ki + . . . + km = n. 

Let Oq be an open neighborhood of the point q G M. We say that a system of coordinates 
ip : Oq ^ W^ is linearly adapted to the flag if, in these coordinates, we have ilj{q) = and 

V-^IA^) =M*^ie...©]R'=% Vi = l,...,m. 

Consider now the splitting R" = M^ © . . .©R ™ and denote its elements x = {xi, . . . , x^) where 
Xi = {x\, . . . , x^') G M'^*. The space of all differential operators in R" with smooth coefficients forms 
an associative algebra with composition of operators as multiplication. The differential operators 



J d 
i = 1, . . . ,m; j = 1, . . . , /cj. We define weights of generators as 



with polynomial coefficients form a subalgebra of this algebra with generators l,x-,^^, where 



d 

1/(1) =0, iyixl) = i, z/(— ^) = -i, 

ox. 



and the weight of monomials 

' dzi ■ ■ ■ dzp 



Ql3 » P 



Notice that a polynomial differential operator homogeneous with respect to v (i.e. whose monomials 
are all of same weight) is homogeneous with respect to dilations 5t : M" — )• M" defined by 

6t{xi,...,Xm) = {txi,t'^X2,...,t"'xm), t>0. (14) 
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In particular for a homogeneous vector field X of weight h it holds 6t*X = t ^X. A smooth vector 
field X G Vec(M"), as a first order differential operator, can be written as 

and considering its Taylor expansion at the origin we can write the formal expansion 

oo 

x^ Yl ^^"^ 

h=—m 

where X^ ' is the homogeneous part of degree h oi X (notice that every monomial of a first order 
differential operator has weight not smaller than —m). Define the filtration of Vec(M") 

pW ={X € Vec(M") : X^^ = 0,Vi < /i}, £ € Z. 

Definition 21. A system of coordinates ip : Oq ^ M" defined near the point q is said privileged for 
a sub-Riemannian structure S if these coordinates are linearly adapted to the flag and such that 
ip*Xi G D'~^) for every i = 1, . . . ,k. 

Theorem 22. Privileged coordinates always exists. Moreover there exist ci,C2 > such that in 
these coordinates, for all e > small enough, we have 

c\ Box(e) C -B(g, e) C C2 Box(e), (15) 

where Box(e) = {x G M", |xj| < e*}. 

Existence of privileged coordinates is proved in P, 1^, O, |l3| . In the regular case the construction 
of privileged coordinates was also done in the context of hypoelliptic operators (see [^). The second 
statement is known as Ball-Box theorem and a proof can be found in |T^. Notice however that 
privileged coordinates are not unique. 

Definition 23. Let S = (M, A,g) be a regular sub-Riemannian manifold and (Xi,...,Xfc) a 
local orthonormal frame near a point q. Fixed a system of privileged coordinates, we deflne the 
nilpotent approximation of S near q, denoted by Sg, the sub-Riemannian structure on R" having 
{Xi, . . . ,Xk) as an orthonormal frame, where Xi := (^*Xj)'~^^ 

Remark 24. It is well known that under the regularity hypothesis, S^ is naturally endowed with 
a Lie group structure whose Lie algebra is generated by left-invariant vector flelds Xi , . . . , X^ ■ 
Moreover the sub-Riemannian distance d in Sg is homogeneous with respect to dilations 6t, i.e. 
d{St{x),5t{y)) = td{x,y). In particular, if Bq[r) denotes the ball of radius r in S^, this implies 

5t{Bq{l)) = Bq{t). 

Theorem 25. The nilpotent approximation Sg of a sub-Riemannian structure S near a point q is 
the metric tangent space to M at point q in the sense of Gromov, that means 

5i/,B{q,e)^Bg, (16) 

where Bq denotes the sub-Riemannian unit hall of the nilpotent approximation Sq. 
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Remark 26. Convergence of sets in ([T6| ) is intended in the Gromov-Hausdorflt topology [12, |2^. In 



the regular case this theorem was proved by Mitchell in |3j]. A proof in the general case can be 
found in lllz 



Definition 27. Let S = (M, A, g) be a regular sub-Riemannian manifold and Sg its nilpotent 
approximation near q. If /i a smooth volume on M, associated to the smooth non- vanishing n-form 
ujfj^, we define the induced volume fiq at the point q as the left-invariant volume on S^ canonically 
associated with u}^{q) G A"'{T*M) (cf. isomorphism (||)). 

From Theorem ^ and the relation^ fi{6£A) = e^fiq{A) + o(e'^) when e — )■ 0, one gets 

Corollary 28. Let ^l he a smooth volum,e on M and Jlq the induced volume on the nilpotent 
approximation at point q. Then, for e — )■ 0, we have 

l^{B{q,e))=e%q{Bq) + o{e'^). 

3.2 Normal forms for nilpotent approximation in dimension < 5 

In this section we provide normal forms for the nilpotent approximation of regular sub-Riemannian 
structures in dimension less or equal than 5. One can easily shows that in this case the only 
possibilities for growth vectors are: 

- dim(M) = 3: g(S) = (2,3), 

- dim(M) = 4: ^(S) = (2, 3, 4) or ^(S) = (3, 4), 

- dim(A/) = 5: ^(S) G {(2, 3, 4, 5), (2, 3,5), (3, 5), (3,4, 5), (4, 5)}. 

We have the following. 

Theorem 29. Let S = (M, A, g) be a regular sub-Riemannian manifold and Sq its nilpotent ap- 
proximation near q. Up to a change of coordinates and a rotation of the orthonormal frame we 
have the following expression for the orthonormal frame of Sq : 

Case n = 3 • ^(S) = (2,3). (Heisenberg) 

Xi = di, 

X2 = 82 + Xi53. 

Case n = 4 • ^(S) = (2, 3, 4). (Engel) 

-^^2 = ^2+ Xid^ + X1X294- 

• ^(S) = (3, 4). (Quasi-Heisenberg) 

X2 = 82 + xidi, 
X3 = 83. 



^Notice that this formula is meaningful in privileged coordinates near q. 
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Case n = 5 • ^(S) = (2, 3, 5). (Cartan) 



X2 = (92 + xid'i + -xld4, + X1X255. 



a(S) = (2,3,4,5). (Goursat rank 2) 

Xi=di, 



^2 = ^2 + xi^s + -xld4 + -x\d^. 
I o 



g(S) = (3,5). ^Coranfc 2) 



Xi = 9i - -X254, 

-^^2 = ^2 + -Xi54 - -Xa^s, 

^3 = ^3 + -X2«95- 



a(S) = (3,4,5). (Goursat rank 3) 



Xl = 9i - -X2di - -X1X295, 
^2 = ^2 + -xidi + -X^(95, 



^3 = 93. 
g{S) = (4,5). (Bi-Heisenberg) 



^1 = 5i - -a;295, 

^2=^2 + -Xi55, 

^3 = 93 - |x455, a G M, (17) 

^4 = 94 + -X355. 

Proof. It is sufficient to find, for every such a structure, a basis of the Lie algebra such that the 
structural constants^ are uniquely determined by the sub-Riemannian structure. We give a sketch 
of the proof for the (2, 3, 4, 5) and (3, 4, 5) and (4, 5) cases. The other cases can be treated in a 
similar way. 

(i). Let S = (G, A,g) be a nilpotent (3,4,5) sub-Riemannian structure. Since we deal with a 
left-invariant sub-Riemannian structure, we can identify the distribution A with its value at the 



*Let Xi, . . . ,Xt be a basis of a Lie algebra g. The coefficients cfj tliat satisfy [Xi,Xj] = '}2i'^ij^i ^^e called 
structural constant of 0. 
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identity of the group Aid- Let {ei, 62, 63} be a basis for Aj^, as a vector subspace of the Lie algebra. 
By our assumption on the growth vector we know that 

dim span{[ei,e2],[ei,e3],[e2,e3]}/Ai(i = 1. (18) 

In other words, we can consider the skew-simmetric mapping 

$(•, •) := [•, •l/A.d : A,d x A^^ ^ T,dG/A,d, (19) 

and condition ([iq ) implies that there exists a one dimensional subspace in the kernel of this map. 
Let X3 be a normalized vector in the kernel and consider its orthogonal subspace D C A^^ with 
respect to the Euclidean product on Aj^^. Fix an arbitrary orthonormal basis {Xi, X2} of D and set 
X4 := [Xi, X2]. It is easy to see that X4 does not change if we rotate the base {Xi, X2} and there 
exists a choice of this frame, denoted {Xi,X2}, such that [X2,X4] = 0. Then set X5 := [Xi,X4]. 
Therefore we found a canonical basis for the Lie algebra that satisfies the following commutator 
relations: 

[Xi,X2] = X4, [Xi,X4] = X5, 

and all other commutators vanish. A standard application of the Campbell-Hausdorff formula gives 
the coordinate expression above. 

(a). Let us assume now that S is a nilpotent (2,3,4,5) sub-Riemannian structure. As before 
we identify the distribution A with its value at the identity and consider any orthonormal basis 
{61,62} for the 2-dimensional subspace Aj^. By our assumption on ^(S) 

dim span{6i,62, [61,62]} = 3 
dim span{6i,62, [61,62], [61, [61,62]], [62, [61,62]]} = 4. (20) 

As in (z), it is easy to see that there exists a choice of the orthonormal basis on Aid, which we 
denote {Xi,X2}, such that [X2,[Xi,X2]] = 0. From this property and the Jacobi identity it 
fojlow^ [X2,[Xi,[Xi,X2]]] = 0. Then we set X3 := [Xi,X2], X4 = [Xi,[Xi,X2]] and X5 := 
[Xi, [Xi, [Xi,X2]]]- It is easily seen that (|20|) implies that these vectors are linearly independent 
and give a canonical basis for the Lie algebra, with the only nontrivial commutator relations: 

[^1,^2] = X3, [^1,^3] = X4, [Xi, X4] = X^. 

{in). In the case (4,5) since dim Ti^G/Aid = 1, the map (|l9| ) is represented by a single 4x4 
skew-simmetric matrix L. By skew-symmetricity its eigenvalues are purely imaginary ±i6i,±i62, 
one of which is different from zero. Assuming 61 7^ we have that a = fo2/&i- Notice that the 



structure is contact if and only if a 7^ (see also Section 5^ for more details on the normal 



form) . n 



Remark 30. Notice that, in the statement of Theorem ^ in all other cases the nilpotent approx- 
imation does not depend on any parameter, except for the (4, 5) case. As a consequence, up to 
dimension 5, the sub-Riemannian structure induced on the tangent space, and hence the Popp 
measure V, does not depend on the point, except for the (4,5) case. 

In the (4, 5) case we have the following expression for the Popp's measure 

where bi, 62 are the eigenvalues of the skew-simmetric matrix that represent the Lie bracket map. 
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Since the normal forms in Theorem 2£ do not depend on the point, except when ^(S) 7^ (4,5), 
we have the foUowing coroUary 

Corollary 31. Let S = (Af, A,g) be a regular sub-Riemannian manifold such that dim(M) < 5 
and G{S) 7^ (4,5). Then if qi,q2 £ M we have that Sg^ is isometric to Sgj o,s sub-Riemannian 
manifolds. 

4 Proof of Theorem [T|: the density is the volume of nilpotent balls 

In this section we prove Theorem ||, i.e. 

UsiQ) = ^MB,). (21) 

It is well known that, being fi absolutely continuous with respect to S^ (see Theorem 20) 



the Radon-Nikodym derivative of /i with respect to S'^, namely f^, can be computed almost 
everywhere as 

KBiq,r)) 



lim „^, „, 
r^oSQ{B{q,r)) 



By Corollary |2^ we get 

fi{B{q,r)) _rQjlg{Bg)+o{rQ) _fig{Bg) 2'^rQ , o{rQ) 



+ 



SQ{B{q,r)) SQ{B{q,r)) 2Q SQ{B{q,r)) SQ{B{q,r))' 

Then we are left to prove the following 
Lemma 32. Let A be an open subset of M. For S^-a.e. q ^ A we have 

l,^S'>iApB(,.r))^^ (22) 

Proof. In the following proof we make use of Vitali covering lemma R and we always assume that 



balls of our covering are small enough to satisfy property (ii) of Proposition 15. 

We prove that the set where (|2^ ) exists and is different from 1 has 5'^-null measure. 
(i). First we show 

S^{Es) = 0, V0<(5<1, 

where 

Es ■.= {qeA:S'^{AnB{q,r)) < (1 - 5)(2r)^,V0 < r < 5} 

Let {Bi} a ball covering of Eg with diam(i?j) < 6 and such that 

Y, dmm{Bif < S^{Es) + e < S^{Es) + e. 



■^Theorem. (Vitali covering lemma, |22[ yj) Let _E be a metric space, B C E and a > such that 'H°'{B) < 00, 
and let J^ be a fine covering of B. Then there exist a countable disjoint subfamily {Vi} C J- such that 

H"{B\[jV.) = 0. 

We recall that J^ is a fine covering of B if for every x £ B and e > there exists V £ J- such that x £ V and 
diam(y) < e. 
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Then we have 

S^{Es)<S'^{An[jBi) 

<^cS«(Ani?,) 

<{l-6){S'^{Es) + e). 

then e — ^ and 1 — 6 < 1 imphes S'^{Es) = 0. 
(ii). Next we prove that 

SQ{Et) = 0, Vt>l, 

where 

Et:={qeA: S'^I.An B{q,r)) > t{2r)^ ,y r small enough}. 

Now let U be an open set such that Et C U and S'^{A nU) < S^{Et) + e. We define 

T := {B{q,r) : q G Et,B{q,r) C U,diauiB{q,r) < 5} . 

Now we can apply Vitali covering lemma to J- and get a family {-Bj} of disjoint balls such that 
S'^iEt \ \Ji B,) = 0. Then we get 

S^{Et)+e>S^{AnU) 

ys'^iAnlJBi) 

>i^diam(S,)^ 

>tS^{Etn[JB^) 

>tSf{Et). 

Letting e, J — )• we have an absurd because t > 1. D 



Since A is open, from this lemma follows formula (|2lD. 

Remark 33. Notice that, for a n-dimensional Riemannian manifold, the tangent spaces at different 
points are isometric. As a consequence the Riemannian volume of the unit ball in the tangent 
space is constant and one can show that it is Cn = vr2/r(^ + 1). Formula @j, where /x = Vol is 
the Riemannian volume, implies the well-known relation between Vol and the (spherical) Hausdorff 
measure 

dVol = ^ d5" = ^ dJT. 

4.1 Proof of Corollary ^ continuity of the density 

In this section we prove Corollary ^ More precisely we study the continuity of the map 

ft.s-q^ MBg)- (23) 

To this purpose, it is sufficient to study the regularity under the hypothesis that p,q does not depend 
on the point. Indeed it is easily seen that the smooth measure /i, which is defined on the manifold, 
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induces on the nilpotent approximations a smooth family of measures {fIq}q£M- In the case fi = V, 
this is a consequence of equahty (^ . In other words we can identify all tangent spaces in coordinates 
with R" and fix a measure /i on it. 

We are then reduced to study the regularity of the volume of the unit ball of a smooth family of 
nilpotent structures in M", with respect to a fixed smooth measure. Notice that this family depends 
on an n-dimensional parameter. 

To sum up, we have left to study the regularity of the map 

q ^ C{%), q G M, (24) 

where Bq is the unit ball of a family of nilpotent structures Sg in M" and C is the standard Lebesgue 
measure. 

Let us denote dq the sub-Riemannian distance in Sg and pq := dq{0, •). Following this notation 



Bq = {x £ R"| Pq{x) < 1} and the coordinate expression (14) implies that 



jCi6aiBq)) = a^CiBq), Va > 0. (25) 

Notice that, since our sub-Riemannian structure is regular, we can choose privileged coordinates 
Tpq : Oq ^ M" smoothly with respect to q. Let now q' / q, there exists a = a{q, q') such that (see 



Remark ^ 



6^Bq> CBqC 6aBq,. (26) 

a 

Using (pSf) , (|26D and monotonicity of the volume we get 



\aQ 7 



C{Bq,) < C{Bq) - C{Bq,) < {a^^ - l)C{Bq 



Then it is sufficient to show that a{q, q') — )■ 1 when q' — )■ q. This property follows from the next 

Lemma 34. The family of functions Pqlx "is equicontinuous for every compact K C M". Moreover 
Pq' ~^ Pq uniformly on compacts inW^, as q' ^ q. 

In the case in which {pt}t>o is the approximating family of the nilpotent distance p, this result 
is proved in [^. See also [^ for a more detailed proof, using cronological calculus. With the same 
arguments one can extend this result to any smooth family of regular sub-Riemannian structures. 
The key point is that we can construct a basis for the tangent space to the structure with bracket 
polynomials of the orthonormal frame where the structure of the brackets does not depend on the 
parameter. 

5 Proof of Theorem |5|: differentiability of the density in the corank 
1 case 

In this section we prove Theorem ^, We start by studying the contact case. Then we complete our 
analysis by reducing the quasi-contact case and the general case to the contact one. 
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5.1 Normal form of the nilpotent contact case 



Consider a 2-step nilpotent sub-Riemannian manifold in M" of rank k. 
Select a basis {^i, . . . , Xf^, Zi, . . . , Zn-k} such that 

A = span{Xi,...,Xfc}, 

n—k 

[Xi,Xj] = ^b'ljZh, i,j = l,...,k, where b'^^ 



V 



-h^. 
"]ii 



(27) 



h=l 



jXi,Zj] = [Zj,Zh] = 0, i = l,...,k, j,h = l,...,n- k. 
Hence the Lie bracket can be considered as a map 

[•,•]: Ax A — >TM/A (28) 

and is represented by the n — k skew-simmetric matrices L = {b^j), h = 1, . . . ,n — k. 

In the contact case we have (fc, n) = {2i, 2i + 1) and our structure is represented by one non 
degenerate skew-symmetric matrix L. Take coordinates in such a way that L is normalized in the 
following block-diagonal form 



L 



fBi 



\ 



\ 



BiJ 



where i?,- := 



-k 
k 



with eigenvalues zizibi,... ,Hbi. Hence we can find vector fields {Xi 
that relations (27) reads 



bi>0. 
,Xi,Yi,... ,Yi,Z} such 



(29) 



A = span{Xi, . . . ,Xe,Yi, . . . ,Ye}, 

[Xi,Y^ = -biZ, i = !,...,£ 

[X^,YJ]=0, i^j 

^[X,,Z] = [Yi,Z] = 0, i = l,...,i 

In the following we call bi, . . . ,b£ frequences of the contact structure. 

We can recover the product on the group by the Campbell-Hausdorff formula. If we denote 
points q = (x, y, z), where 

x = {xi,...,Xi)£R\ y = {yi,...,ye) £R\ zeR, 

we can write the group law in coordinates 



q-q = [x + x'^y + y'^z + z'--^ bi{xix[ 



yiVi) 



(30) 



Finally, from (|30|), we get the coordinate expression of the left-invariant vector fields of the Lie 
algebra, namely 

Xi = dx^ + -biyidz, i = !,...,£, 

li Oy^ '^'^iXiOzi 1 J-) • • • )^) 

z = d,. 

In this expression one of frequences bi can be normalized to 1. 



(31) 
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5.2 Exponential map in the nilpotent contact case 

Now we apply the PMP to find the exponential map Exp where qq is the origin. Define the 
hamiltonians (linear on fibers) 

K^{X) = {X,X,{q)), h,,{X) = {X,Yi{q)) , K{X) = {X,Z{q)). 

Recall from Section [2.2| that q{t) is a normal extremal if and only if there exists \{t) such that 

Xi — tli 

Vi = Vi (32) 

where 

u,{t) := /i„,(A(t)), v,{t) := K^{X{t)), w{t) := /i^(A(t)). 

Remark 35. Notice that from ( |32D it follows that the sub-Rieniannian length of a geodesic coincide 
with the Euclidean length of its projection on the horizontal subspace {xi, . . . ,Xn,yi, ■ ■ ■ , Vn)- 




m = f {jZ^ui{t)+vi{t))\\t. 



Now we solve ( ^2| ) with initial conditions (see also Remark 13) 

(xO,y°,zO) = (0,0,0), (33) 

(w°, 7;°, w^) = {ul, ...,4,vl,..., v'l, w^) G S^^-^ X M. (34) 

Notice that w = w^ \s constant on geodesies. We consider separately the two cases: 
(i) If ui / 0, we have (denoting ai := hiw) 

Ui{t) = Uj cos ait — Vi sin ait, 

Vi{t) = Ui sinajt + Vi cos ait, (35) 

w(i) = w. 



Prom (32) one easily get 



Xi (t) = — (n^ sin ait + v^ cos a^t — v^), 
ai 

yi{t) = — {—u^ cos ait + v^ sin ait + u^) , (36) 

a,; 



W = ^(^* - E r((^°)' + K°)')) si^«^*)- 



2u;2 
{ii) If u; = 0, we find equations of straight lines on the horizontal plane in direction of the vector 

Xi{t) = u'^t yi{t) = v^t z{t) = 0. 
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Remark 36. To recover symmetry properties of the exponential map it is useful to rewrite ( pq ) in 
polar coordinates, using the following change of variables 

Ui = riCos6i, Vi = risin9i, i = !,...,£. (37) 



In these new coordinates (36) becomes 



r. 



Xi{t) = — (cos(ojt + 9i) — cos 6i), 
ai 

yi{t) = — (sin(ajt + 6i) - sin^j), (38) 

Cli 



1 yl' 



2u;2 

t 

and the condition (ii°, v^) G S"^^'^ implies that r = (ri, . . . , r^) e S*^. 

Prom equations ( |38| ) we easily see that the projection of a geodesic on every 2-plane {xi,yi) is 
a circle, with period Tj, radius pi and center Cj where 

Ti = - — , Pi = —^ d = ---^{cos6i,s'm6i), Vi = !,...,£ (39) 

biW biW biW 

Moreover ( |32|) shows that the z component of the geodesic at time t is the weighted sum (with 
coefficients bi) of the areas spanned by the vectors (xj(t),yj(t)) in M^. 

Lemma 37. Let 7(i) be a geodesic starting from the origin and corresponding to the parameters 
{ri,0i,w). The cut time tc along 7 is equal to first conjugate time and satisfies 

■ '' (40) 



w maxj bi ' 

with the understanding that tc = +00, if w = 0. 

Proof. We divide the proof into two steps. Recall that a geodesic lose optimality either at a cut 
time or at a conjugate time, see Definition 17. First we prove that (HO) is a conjugate time and 



then that for every t < t^ our geodesic is optimal. 

The case if = is trivial. Indeed the geodesic is a straight line and, by Remark ^, we have 
neither cut nor conjugate points. Then it is not restrictive to assume w ^ 0. Moreover, up to 
relabeling indices, we can also assume that bi > b2 > ■ ■ ■ > bi > 0. 

(i) Since, by assumption, 5i = maxj6j, from (|38|) it is easily seen that projection on the (xi,yi)- 
plane satisfies 

xi{tc) = yi{tc) = 0. 

Next consider the one parametric family of geodesic with initial condition 

(ri,r2,...,r£,6'i + (/>, 6'2, . . . , 6'^,'u;), cf) G [0,27r]. 

It is easily seen from equation (|3g) that all these curves have the same endpoints. Indeed neither 
{xi,yi), for i > 1, nor z depends on this variable. Then it follows that tc is a critical time for 
exponential map, hence a conjugate time. 
{ii) Since w ^ 0, our geodesic is non horizontal (i.e. z(t) ^ 0). By symmetry, we can focus on the 
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case w > 0. We know that, for every i, the projection of every non horizontal geodesic on on the 
plane {xi,yi) is a circle. Moreover for the z-th projected curve, the distance from the origin is easily 
computed 

biwt, , . , X sinx 



r]i{t) = \/xi{tY + yi{tY = rit sinc(^-— ), where sinc(x 



2 L.\-~/ 

X 

Let now i < tc, we want to show that there are no others geodesies 7(t), starting from the origin, 
that reach optimally the point 7(t) at the same time t. Assume that 7(t) is associated to the 
parameters {ri,0i,w), where (fi, . . . ,r^) E S^, and let us argue by contradiction. If j{t) = j{i) it 
follows that rji^t) = rji{t) for every i, that means 

,biwi _ _ ,biwi 
nt sulci— —) =rit smc{——), ^ = !,...,£. (41) 

Notice that, once w is fixed, fj are uniquely determined by ( |4l| ) (recall that i is fixed). Moreover, 
9i also are uniquely determined by relations (|3^). Finally, from the assumption that 7 also reach 
optimally the point 7(t), it follows that 

27r bjwi , , „ , , , 

t<tc = —= =^ ^<^ Vi = l,...,l 42 

Since sinc(x) is a strictly decreasing function on [0,7r] it follows from (^Tj) that, ii w ^ w we have 
w>w =^ fj>rjVi = l,...,f =^ ^ ^rf > ^ ^r^^ = 1, 

j i 

w < w =^ rj<rjVz = l,...,f =^ >^ fj < N^ r^ = 1, 
which contradicts the fact that (fi, . . . ,?£) G 5^. 



D 



Consider now the exponential map from the origin 



Xi = - — (cos(biW + 6i) — cosOi 
biw 



{r,0,w) I— )■ Exp(r, 6*,^;) = < 



= -^(sin(6iu; + Oi) - smOi) (43) 

biW 

- ^i^\r\' - E. ^ sin6,u;) = E. ^2(b^^ - ^^^b.^)- 

By Lemma ^, the set -D where geodesies are optimal and their length is less or equal to 1, is 
characterized as follows 

D := {(r, 9, w), \r\ < 1, \w\ < 27r/max6j}. 

Thus the restriction of the exponential map to the interior of D gives a regular parametrization of 
the nilpotent unit ball B, and we can compute its volume with the change of variables formula 



C{B)= idC= / \AQiJ\Rdrdedw, (44) 

where J denotes the Jacobian matrix of the exponential map (^). Notice that we have to integrate 
with respect to the measure dudvdw = RdrdOdw, where R = J^j ?^j because of the change of 
variables (37). 
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Lemma 38. The Jacobian of Exponential map ( |43| ) is given by the formula 



det J{r, 6, w) 






where we denote with B = W^ hi . 

Proof. We reorder variables in the following way 

{ri,ei,...,re,0i,w), {xi,yi, . . . ,xe,yi, z), 

in such a way that the Jacobian matrix J of the exponential map ( p3|) is 

/Qi 



J 



Q2 



\Zl Z2 



Wi\ 

W2 

Qi We 
Zi du,zj 



where we denote 



^ ' " \dr,yi de^yij 

Vr,:=(^-^^), Zr.= {dr^z de^z) 

\^wyij 



i = 1 



1 • • • 1^1 



Notice that xi^yi depend only on rj,^j. 

To compute the determinant of J, we write 



z = z{r,w) =^^Zi{ri,w), Zi :-- 



2bjW 



:{biW — sinbiw) 



and we split the last column of J as a sum 



/Wi\ 

W2 




( Wi\ 



+ 


( ' \ 

W2 


Wi 
\dyjz) 






Kd^zJ 






\dwZ2/ 



+ ••• + 



/ \ 





We 
\dwzej 



(45) 



(46) 



(47) 



(48) 



Notice that in the i-th column only the i-th variables appear. 

By multilinearity of determinant, det J is the sum of the determinants of i matrices, obtained 
by replacing each time last column with one of vectors appearing in the sum (|48|). If we replace it, 
for instance, with the first term, we get 



/Qi 



J^ 



Q2 



Wi\ 




Qe 
\Zi Z2 ... Ze du,zi) 
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Now with straightforward computations (notice that dgzi = in Zi), we get 



Setting 
we find 



det Ji = det Q2- ■ ■ det Q^ ■ det , „ 

Zi OwZl 

detQ2---detg,-det( „^'i ^i ^^ ,. 



A^ = {Q'i Wi) i = l,...,l, 



det Ji = det Q2- ■ ■ det Qi ■ {d^zi det Qi + dr^zi det Ai). 
Similarly, we find analogous expressions for J2, ■ ■ ■ , Ji- Then 

e e 

det J = y^ det Ji = /^( I I det Qj){dwZi det Qi + 9,.; 2* det Ai 



From (|43D, by direct computations, it follows 




Moreover, with some computations, one can get 






o ''^f ■ 2M'^\ '^f tu ■ u \ 

OwZi = -^ sm (— -) - 7 — :T[OiW - smbiW) 
w^ 2 biW^ 

^riZi = T^ibiW - sinhiw) 
biW^ 

Arf . biW fbiw Mw . Mw 
d«t^^ = -^2^^^^( — ) ( — ^°^( — ) - ^^"( — ) 



4r? . ,biW ^ f biW ,biW , . ,biW . 



and finally, after some simplifications 

dyjZideiQi + dr^ZideiAi = --p^sin(^) ( ^cos(^) - sin(- 
from which we get (|45|). D 



From the explicit expression of the Jacobian (H3) we see that integration with respect to hori- 
zontal variables {ri,6i) does not involve frequences, providing a costant Ci. Hence, the computation 
of the volume reduces to a one dimensional integral in the vertical variable w: 

fT^ Ce sr^ Tf ■ 2,bjw\ . hw fbiW Mw . Mw \ 

y = y_^^ -^2^-2ir2 L n ^- (— ) «-(— ) [— ^°«(— ) - ^-(— )J ^^- 

maxfe^ 2 = 1 \j^i J 

Using simmetry property of the jacobian with respect to w and making the change of variable 
2s = w, the volume become (reabsorbing all costants in Ci) 

y ^ ' "H2~2i+2 X] JJsin^(6js) 1 sin(6js)(6iScos(6.js) - sin{bis))ds. (49) 
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Remark 39. In the case i = 1 (Heisenberg group) we get 

V = —{1 + 2tt Si(2^)) ~ 0.8258, Si(x) := / dt. 

1^ JO * 

Notice that this is precisely the value of the constant f-ps^ since in this case Popp's measure coincide 
with the Lebesgue measure in our coordinates. 

5.3 Differentiability properties: contact case 

Let us come back to the differentiability of the map 

q^C{Bg). (50) 

A smooth family of sub-Riemannian structures is represented by a smooth familly of skew 
symmetric matrices L{q) (see Section q]^). Recall that, for a smooth family of skew-symmetric 
matrices that depend on a n-dimensional parameter, the eigenvalue functions q i— )• bi{q) exists and 



are Lipschitz continuous with respect to q (see [31 1). 



Thus, if we denote with V{q) the volume of the nilpotent unit ball corresponding to frequences 
bi{q), . . . , ht,{q), formula ( p9|) can be rewritten as 

V{q)= n Giq,s)ds, (51) 

Jo 



where 



a{q) :-- 



max bi (q) ' 
G{q,s) := 2^+2 R2 X] n^™^(^j(^)'^) sm{bi{q)s){bi{q)s cos{bi{q)s) - s'm{bi{q)s)). (52) 






Notice that we have dropped the constant Ci that appear in ([49|) since it does not affect 
differentiability of the volume. 

Remark 40. Since the family of sub-Riemannian structures q i— )■ L(q) is smooth, the exponential 
map smoothly depends on the point q. As a consequence the integrand G{q, s), being the Jacobian 
of the exponential map, is a smooth function of its variables. 

In addition, altough q — t- bi{q) is only Lipschitz, it is easy to see that the function a{q) = 



max6,(g) is semiconvex with respect to q (see also [18|). In particular a{q) admits second derivative 
almost everywhere. 

If for all q all bi{q) are different (i.e. there are no resonance points), then the eigenvalue functions 
can be chosen in a smooth way. As a consequence the volume V{q) is smooth, since all functions 
that appear in (^Tj) are smooth. This argument provides a proof of Proposition ^ in the case of 
contact structures. 

On the other hand, we prove that, along a curve where a{q) is not smooth, i.e. when the two 
bigger frequences cross, V{q) is no longer smooth at that point. 
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Proof of Theorem [^. In this proof by a resonance point we mean a point qq where the two (or more) 
biggest frequences of L(go) coincide. More precisely, if we order the frequences as 

bi{q)>b2{q)>...>be{q), 

a resonance point is defined by the relation bi{qo) = b2{qo)- As we noticed, V is smooth at non- 
resonance points. 

We divide the proof into two steps: first we prove that V (i C^ and then we show that, in 
general, it is not smooth (but C^ on smooth curves). 

(i). We have to show that V is C^ in a neighborhood of every resonance point qQ. First split 

the volume as follows 

ra{<lo) ra(q) 

V{q)= G{q,s)ds+ G{q,s)ds (53) 

JO J a{qo) 

The first term in the sum is smooth with respect to q, since it is the integral of a smooth function 
on a domain of integration that does not depend on q. We are then reduced to the regularity of 
the function 

ra{q) 

W{q):= G{q,s)ds (54) 

Ja{qo) 

We have the following key estimate 



Lemma 41. Let qo £ M be a resonance point. Then, for any neighborhood ofqo, there exists C > 
such that 

ra.{q) 

G{q, s)ds 

'a(qo) 



<G\q-qo\^ (55) 



Proof. It is sufficient to prove that every derivative up to second order of G vanish at ((70) ^^(^o))- 
Indeed, being G smooth, computing its Taylor polynomial at {qo,a{qo)) only terms with order 
greater or equal than three appear (both in g — go and s — a{qo)). Thus, integrating with respect 
to s and using that \a{q) — a{qQ)\ = 0{\q — qo\), we have the desired result. 

From the explicit formula ( p2|) it is easy to see that G{q, a{q)) = for every q G M. In particular 
G{qo, a((7o)) = 0. Moreover, since at a resonance point qo at least the two bigger eigenvalues coincide, 
say bi{qo) = b2{qo) = /3, we have a{qo) = n/fS and in a neighborhood of {q, s) = {qo, a((7o)) 

sin2(6i(g)^) sin(62((?)^) = OiMq) - b^iq)^) = Oi\q - qo\'). (56) 

due to the Lipschitz property of bj{q), for j = 1, 2. 

From (m ) and (pq ) one can easily get that every derivative of G up to second order (in both 
variables q and s) vanish at {qo,a{qQ)). 

a 

To show that V £ C^ we compute the first three derivatives of M^ at a non resonant points q, 
and we show that, when q tend to a resonance point qo, they tends to zero. We then conclude the 
continuity of the first three derivatives by Lemma ^. 

In the following, for simplicity of the notation, we will denote by ^ the partial derivative with 

respect to some coordinate function on M. For instance ^-^ denote some second order derivative 
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At non resonance points q we have 



-{q) = G{q, a{q)) Trig) + / ^^(^^ s)ds 



dq . 'dq Ja(qo) dq 

=0 

-^-{q,s)ds 

a(go) ^9 



which tends to zero for g — >• go- Using Lemma 41 we conclude V £ C^. Next let us compute 

d'^W , . dG, , ..da, . f<'i> d'^G 

dq^ uq uq Ja(qo) 



dG da /■'^(s) g'^Q 



Since ^{qo,a{qo)) = (see proof of Lemma |4l|) and ^ is bounded by Lipschitz continuity of the 



maximum eigenvalue, it follows that ^-^ tends to zero as g — )■ gO) and using again Lemma ^ we 
have that V € C'^. In analogous way one can compute the third derivative 

d^W , , d^G , .^^/9a,^^^ d'^G, , ,,da , , 



dq^ dqds ' \dq J dq^ ' dq 

d'^a , . r^'i^ d^G 



dG da /■'^w) Q-iQ 



i{go) 

dq 



Using again that every second derivative of G vanish at {qo,a{qo)) and that ^ is bounded, it 
remains to check that ^{q,a{q))i-^{q) tends to zero as g — ?• go- From (pq ) one can see that 

^ = 0{\bi{q) - 62(g)P)- Hence it is sufficient to prove that = 0{l/\bi{q) - 62(9)!), which is a 
consequence of the following lemma. 

Notice that, for every skew-symmetric matrix A, with eigenvalues ziziXj, with j = 1, . . . ,n, the 
matrix iA is an Hermitian matrix which has eigenvalues ±Aj, j = 1, . . . , n. 

Lemma 42. Let A, B be two n x n Hermitian matrices and assume that for every t the matrix 
A + tB has a simple eigenvalue Xj{t). Then the the following equation is satisfied 

.■_^^\_{^^ (59) 



k^j 



Xj — Xk 



where {xk{t)}k=i,...,n is an orthonormal basis of eigenvectors and Xj{t) is the eigenvector associated 
to Xj{t). 

Proof. In this proof we endow C" with the standard scalar product {z, w) = X]fc=i z0Jrk. Since 
Xj{t) is simple for every t, both Xj{t) and the associated eigenvector Xj{t) can be choosen smoothly 
with respect to t. By definition 

{A + tB)xj{t) = Xj{t)xj{t), \xj{t)\ = 1. 

Then we compute the derivative with respect to t of both sides 

{A + tB)xj{t) + Bxj{t) = Xj{t)xj{t) + Xj{t)xj{t), (60) 
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and computing the scalar product with Xj{t) we get 

Xj{t) = {Bxj{t),Xj{t)) , hence Xj{t) = 2Re {Bxj{t),Xj{t)) . (61) 

using that A + tB is Hermitian and {xk{t),Xk{t)) = 0. On the other hand, the scalar product of 
with Xk , with k ^ j gives 

[Xj, Xjg) 



Xj — Afc 
Substituting Xj = Ylk=i i^jj^k) Xk in (|6l|) we have (|59|). D 



Then the third derivative goes to zero for q ^ qo, and using again Lemma 41 we conclude that 



(n). Now we study the restriction of the map (pOl ) along any smooth curve on the manifold, and 
we see that, due to a simmetry property, V is C^ on every curve but in general is not C^. Notice 
that (|55|) gives no information on the fourth derivative but the fact that it is bounded. Indeed it 
happens that it is continuous on every smooth curve but its value depend on the curve we choose. 

From now on, we are left to consider a smooth one-parametric family of sub-Riemannian struc- 
ture, i.e. of skew symmetric matrices. 

Remark 43. An analytic family of skew-simmetric matrices 1 1— )• L{t) depending on one parameter, 
can be simultaneously diagonalized (see again [^l|), in the sense that there exists an analytic (with 
respect to t) family of orthogonal changes of coordinates and analytic functions bi{t) > such that 

'0 -biity 



L=\ •.. |, where B.it) := ^^^^^^ -)• (62) 

Be{t)^ 

In the case of a C°° family 1 1— t- L{t), we can apply the previuos result to the Taylor polynomial of 
this family. As a consequence we get an approximate diagonalization for L{t), i.e. for every A^ > 
there exists a smooth family of orthogonal changes of coordinates and smooth functions bi{t) > 
such that every entry out of the diagonal in L{t) is o{t ). Namely 

/B,{t) o{t^)\ 

m=\ •.. 1 where B.it) := (^^^ " ^If ) . (63) 

W) Beit)) ^'" ' '^ 

Since we are interested, in the study the C^ regularity of (^), for k finite, in what follows we 
can ignore higher order terms and assume that L{t) can be diagonalized as in the analytic case 

I- 



Prom the general analysis we know that V is C^. To prove that t i— )• V{t) is actually C^ we 
discuss first the easiest case i = 2 (i.e. the contact (4,5) case) and then generalize to any i. 

(i) Case i = 2. To start, assume that bi{t), b2{t) cross transversally at t = . This means that 
for the volume V{t) we have the expression 



V{t) 



r''' G{t,s)ds, ift>0, 6,(0=62(0) 

•^V and (64) 

Jo 
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Since the regularity of the volume does not depend on the value 5i(0) = 62(0), we can make the 
additional assumption 

bi{t) = l + tci{t), i = l,2 

for some suitable functions ci(i),C2(t). Notice that a'{t) is discontinuous at i = and the left and 
right limits are 

a', := lim a'(t) = — vrcifO), a'_ := lim a'(t) = — vrcofO), 

From the explicit expression of G it is easy to compute that 

-(0, a(0)) = ^ciC2(ci + C2), 



(0, a(0)) = ^(c? + 4ciC2 + 4), ir^iO, a(0)) = ^(ci + C2), 



d'^G ,„ ,„,, 2,2, 2^ ^^G' ,„ ,„,, 6 



dt^ds''''' ^3,^1 , ,.,.z , -2y> g^^^2. ... ^. 

where we denote for simplicity q := Cj(0). 

Let us compute ^^- To this purpose let us differentiate with respect to t formula (|58|) (where 
q is replaced by t). Using the fact that all second derivatives of G vanish at {t,s) = (0, a(0)) (see 
the proof of Lemma ^) we have that the 4-th derivative of VF at t = is computed as follows 



where W is defined in ([5^). It is easily checked that W^^' is continuous (but does not vanish!). 
Indeed we have 

lim W^^\t) = lim W^^Ut) = -—44. 
t-5>o+ t^o- vr 

The same argument produce an example that, in general, V{t) is not C^. Assuming 

bi{t) = 1 + tcj, ci / C2 constant, i = 1,2, 

a longer computation, but similar to the one above, shows that 

lim W'^^^t) = c?(13c? - 29ciC2 + 22c^) 

i->-0+ IT 

lim W^^Ut) = ci(13c^ - 29ciC2 + 22c?) 

t-s-O- vr 

and the 5-th derivatives do not coincide. 

Remark 44. Notice that the assumption of transversality on bi{t) and 62 (^) at t = is not restrictive. 
Indeed if 61 (t) — &2(0 = 0{t ) for some A; > 1, then from the proof of Lemma ^ (see in particular 
(|56|)) it follows that at least 2k + 1 derivatives of G vanish at (t, s) = (0, a(0)), increasing the 
regularity of W. 
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{ii) General case. We reduce to case (i). 

We can write G(t, s) = X]j=i Gi{t, s) and V(t) = X]i=i ^(^) where we set 



G^{t,s) 



s2£+2 



1 rsin^(6j(i)s) sm{bi{t)s){bi{t)scos{bi{t)s) — sm{bi{t)s)), 

a{t) 



(65) 



Vi{t):= I G^{t,s)ds, i = l,...J. 

Jo 

Assume that bi, 62 are the bigger frequences and that they cross at t = 0, i.e. 

bi{t) <b2{t) <bi{t), Vt<0, Vi = 3, ...,n. 
bi{t) <bi{t) <b2it), Vt>0, Vz = 3, ...,n. 

Prom the exphcit expression above it is easy to recognise that for Gi and G2 we can repeat the 
same argument used in (i). Indeed if we denote with G{t, s) the integrand of the (4, 5) case we can 
write Gi + G2 as the product of a smooth function and G 

Gi{t, s) + G2{t, s) = I ^ n sin2(6,(i)s) j Git, s), 

which imphes that that Vi + V2 is a C^ function. 

Moreover it is also easy to see that V3, . . . ,Vn are C^. Indeed from the fact that 61 (t) and 62(0 
both appear in sin^ terms in Gi{t, s) for i > 3, it follows that in this case 

f)G rfiC- r)^G 

G,{t,a{t))^^{t,a{t))^0, ^(0,a(0))=0, ^(0,a(0))=0, i = 3,...,n, 



dt 



df^ 



dt^ 



and we can apply the same argument used in (i) to the function ¥({€) = j^^^' ^^{t,s)ds, for 



r«(t) dGi , 



,n, showing that it is C , that means Vi ^ C 
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Remark 45. As we said the value of the 4-th derivative depend on the curve we chose, hence we 
cannot conclude that V £ C^ in general. Moreover, we explicitly proved that V ^ C°° since in 
general is not C^, even when restricted on smooth curves. 

Moreover from the proof it also follows that, if more than two frequences coincide at some point 
(for instance if we get a triple eigenvalue) , we have a higher order regularity for every Vi , and the 
regularity of V increases. 

5.4 Extension to the quasi-contact case 

Recall that in the quasi contact case the dimension of the distribution is odd and the kernel of the 
contact form is one dimensional. Hence, applying the same argument used in Section 5^, we can 
always normalize the matrix L in the following form: 



L 



fB, 



V 



Bf 



\ 



0/ 



B, 



b, > 0. 



30 



In other words we can select a basis {Xi, . . . , Xi, Yi, . . . , Yi, K, Z} such that 

'A = spanjXi, . . . ,Xi,Yi, . . . ,Yi, K}, 

[X„Y^ = -b,Z, i = !,...,£ 

[Xi,Yj]=0, i^j (66) 

[Xi,K] = [Yi,K] = 0, i = l,...J 

[[Xi,Z] = [Y,,Z] = 0, i = l,...,i 

where the new vector field K is in the kernel of the bracket mapping, i.e. it commutes with all 



others elements. Since abnormal extremals are never optimal in quasi contact case (see Remark 14), 
we are reduced to compute the exponential map to find geodesies. With analogous computations 
of contact case we get the following expression for the exponential map from the origin 



Xi{t) 

Viit) 



— ^(cos(6jtt;t + 6i) — cos 6i) 
biw 



biw 

X2£+lit) =U2e+lt, 



{sm{biwt + 6i) — sinOi), 



(67) 



z{t) 



1 "T^r"^ 

— {\r\'^wt-}2irsmb^wt). 



2w'^"' ' "" ^^ bi 

i 

Prom (|67|) it is easily seen that the jacobian of exponential map has exactly the same expression 



as in contact case (43). Since zero is always an eigenvalue of L, but is never the maximum one, we 
can proceed as in the contact case and all the regularity results extend to this case. 



6 Proof of Theorem g and extension to general corank 1 case 

We start this section with the proof of Theorem |^, after that we extend the result to the general 
corank 1 case. 

Proof of Theorem^. Let S be a sub-Riemannian structure such that dimAf < 5. 

(i). If ^(S) 7^ (4,5). From Theorem 29 we know that at every point q G M, the nilpotent 
approximation Sg has a unique normal form, hence by Corollary ^ all nilpotent approximations 
are isometric. From this property it easily follows that /-ps, the Popp volume of the unit ball, is 
constant (recall that Popp measure is intrinsic for the sub-Riemannian structure). This also implies 
that for a smooth volume /x the density f^s is smooth. 

(ii). If ^(S) = (4,5) by Theorem |2^ it is sufficient to consider the case when the family of 
nilpotent structure has the normal form (|l7|), where a = a{q) depends on the point. Notice that 
the formula (4' 

Theorem 



for the volume of the unit ball is still valid, where now 6i = 1 and 62 = l^l- 
5| proves that the density is C'^ at points where \a\ > 0, i.e. in the contact case. We 
are then reduced to the study of the volume near a point where the eigenvalue a crosses zero. In 
particular we show that the volume is smooth at these points. Since the eigenvalue a is approaching 
zero, it is not restrictive to assume \a\ < 1. Let us consider then the function defined on the interval 
(-1,1) 



W{a) 



1 



a2s6 



(sin (as) sin s(s cos s — sins) + sin ssin(Q;s)(ascos(as) — sin(Qs))) ds. (68) 
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Note that V{a) = W{\a\), where V{a) denotes the volume of the nilpotent bah relative to frequences 
1 and a < 1. It is easy to see that both 



sin^ as , 1 



and — ^ sin(as)(ascos(as) — sin(as)) (69) 



a^ a"^ 



are smooth as functions of a (also at a = 0). Hence P^ is a smooth function (for a G (—1,1)). 
Moreover it is easy to see that W is an even smooth function of a. Thus W it is smooth also as a 
function of |a|, which completes the proof. 

The same argument applies to prove that the C^ regularity holds in the general corank 1 case. 



Indeed, from (52) and the fact that (|69| ) are smooth functions at a = 0, it follows that the integrand 
G{q,s) is smooth as soon as one of the eigenvalue bi(q) is different from zero (recall that bi > 
by definition). Since the structure is regular (i.e. the dimension of the flag do not depend on the 
point) and bracket generating, we have that maxj bi{q) > for every q, hence the conclusion. D 
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